Macroscopic quantum phases of a deconfined QCD matter at 

finite density 



Jin Hosek 

Institute of Nuclear Physics, Rez near Prague, CZ 250 68, Czech Republic 

(February 1, 2008) 



00 

<3\ 



Abstract 



O . Formalism for a unified description of distinct superfluid phases of a decon- 

J-n \ fined QCD matter at finite density together with the phase of spontaneously 

broken chiral symmetry is presented. Dispersion laws of the quasiquark ex- 
' citations in both diamagnetic and ferromagnetic phases with spontaneously 

broken chiral symmetry are exhibited explicitly 

> 

m 

At large enough baryon number densities the hadrons overlap. Due to the asymptotic 
freedom of QCD [l[ in which we trust such a hadronic system should behave as a weakly 
interacting quark-gluon plasma [[J. Restricting ourselves to temperatures at which the 
thermal wavelength becomes comparable with interparticle spacing we can expect behavior 
Oh! of the system characteristic of a genuine quantum Fermi liquid. 

Ph| In this paper we consider the chiral SU(2) limit of QCD which is a good approximation 

to the real world within a 20% accuracy. Numerical value of the quark number density 
n = \pp at T = at which the deconfinement sets in can be estimated by taking n = 
0.72/m~ 3 , the quark number density inside the nucleon. Coherence temperature at which 
the quantum behavior sets in is estimated as follows. Average energy p of a massless quark 
equals n c rifn s 3^k B T which defines the thermal wave length X T = \^f- It should be larger 
than n _1//3 , the average interquark spacing. Hence, below the line T = ^Tm 1 ^ 3 (with h = kg 
set equal to one) in the (n^ 3 ,T) QCD phase dia gram we expect the behavior of the system 
characteristic of a genuine massless Landau Fermi liquid. For example, its specific heat 

Cy ~ T. 

If at the Fermi surface of a Landau Fermi liquid there is an attraction between two 
fermions with momenta pp and— pp respectively, the system becomes unstable with respect 
to the spontaneous condensation of Cooper pairs. This happens in the electron Fermi liquids 
in metals, and turns that system into a BCS superconductor [[|. This happens also in the 
Fermi liquid of 3 He, and turns that system into a chain of distinct superfluid phases M. It 
is clearly important to ask : what are the realistic interactions in the QCD Fermi liquid, 
whether it is of the Landau type, and whether the Cooper instability takes place also there. 

Unfortunately, there are no experimental data at present, either real or the lattice ones 
which would check our considerations. There are, however, rather solid theoretical arguments 
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that to a reasonable approximation the effective interactions between massless quarks in the 
Landau QCD Fermi liquid are governed by an SU (3) c x SU (2) L x SU (2) R globally invariant 
local four-quark Lagrangian Ci nt . The individual terms in £ int originate from : (i) massive 
(due to Debye screening) chromoelectric gluon exchange; (ii) massive (due to dynamical 
Higgs mechanism) chromomagnetic gluon exchange; (iii) quark interactions with instanton 
zero modes; (iv) dynamically generated massive collective excitation exchange. Since at 
least points (ii) and (iv) are not under a safe theoretical control it is necessary to have a 
method capable of analyzing the phase structure of the Landau QCD Fermi liquid for a 
general form of C int . 

Origin of the Cooper-pair instability is understood at present within the effective field- 
theory framework both in nonrelativistic f| and relativistic [U interacting Fermi systems. 
The local four-fermion couplings listed above do contain the attractive channels necessary for 
Cooper pairing, and the renormalization group determines their running from the matching 
point to the infrared. This provides another good reason for developing a general method 
of analysis of the phase structure of a low-temperature system of many fermions carrying 
spin, flavor and color. With perturbative forces in mind the idea of color superfluidity was 
mentioned already in 0, and elaborated in 0. Recent papers PHlOl deal already with local 
four-quark interactions discussed above. 

With a general local four-fermion interaction between quarks carrying spin, isospin and 
flavor there can be four types of the superfluid ground-state condensates different from zero 
due to Pauli principle (here we ignore the alluring possibility of spontaneous parity violation 
within parity-conserving QCD) : 

V(i) = (^Ar 275 ^ c ) , (1) 

V(?) = (?Mt 3 t 2 7o73^ C ) ' ( 2 ) 
V (3) = (^<SV 3 t 2 75^ C ) , (3) 

f(4) = (^^27073^) (4) 

Although we use the convenient Lorentz-covariant notation with ip c = Cip where C is the 
charge-conjugation matrix it is clear that the only sacred property of the ground state here 
is the translation invariance. We contemplate four distinct macroscopic quantum phases : 

1. Condensate vn\ with the color-antisymmetric Clebsch-Gordan (CG) matrix A chosen 
as A ab = ie ab3 = — (A 2 ) ab corresponds to a ground-state expectation value of the order 
parameter <3> c . It describes a Lorentz scalar, isosinglet color triplet superfluid. 

2. Condensate v^) corresponds to a ground-state expectation value of the order parameter 
^ c j. fiu - It describes a color-triplet superfluid which at the same time behaves as ordinary 
as well as flavor ferromagnet : Qj.^ is an isospin 1 Lorentz antisymmetric tensor field. 

3. Condensate t>(3) with the color-symmetric CG matrix S chosen as S ab = |5 ab — ^(Ag) ab 
corresponds to a ground-state expectation value of the order parameter It de- 
scribes a Lorentz scalar color-sextet superfluid which at the same time behaves as a 
flavor ferromagnet. 
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4. Condensate t>( 4 ) corresponds to a ground-state expectation value of the order parameter 
It describes an isoscalar color-sextet superfluid which at the same time behaves 
as an ordinary ferromagnet. 

Recent papers PHTtl are devoted predominantly to the first case with some estimates || 
made also of the fourth possibility. 

It is important to realize that all terms (|1}^) are invariant with respect to global chi- 
ral rotations. This implies that they cannot produce physical, chiral-symmetry-violating 
quark masses in the deconfined phase. There is, however, a ground-state quark- ant iquark 
condensate 

w = (ipip) (5) 

which can do namely this, and it should be considered together with the quark-quark con- 
densates It is well established that at finite density the local four-quark interactions 
considered above can give rise to the quark masses IPIBI and, by virtue of the Goldstone 
theorem also to the massless pions. 

In this Letter we present main steps of a generalization of the field-theory approach 
developed for superconductivity by Nambu |T^], and for spontaneous breakdown of chiral 



symmetry in the nucleon-pion system by Nambu and Jona-Lasinio ||13|| . It gives almost 
immediately the dispersion laws of the true fermionic excitations as exemplified explicitly 
below. It enables the (numerical) analysis of the phase structure of the system for a given 
set of couplings, as will be clearly seen from the following. It also provides a systematic 
way of investigating the gapless collective Nambu- Goldstone excitations corresponding to 
all symmetries spontaneously broken by the condensates This part of the program 

will be published separately. 

A self-consistent perturbation theory is defined by its bilinear Lagrangian with all terms 
(|1|-||) taken into account : 

C' = C - £ s - C A = $ + /i7o)^ - ipZip - 

-\tP(Aj 5 t 2 \i) + A7o73r 3 T2A ( 2) + Sj 5 t 3 t 2 \ 3 ) + S^ 3 T 2 A {i) )i) C + H.c.) (6) 

It is assumed that its ground state for Auj, £ different from zero is energetically advantageous 
with respect to the naive ground state corresponding to £ = $ + Hlo)^- Nonzero values 
of A(j) and E are eventually found by solving the gap equations expressing the condition that 
the lowest-order perturbative contribution of C int = C int + £a + using the propagator 
defined by £ , vanishes. 

Main trick which enables to realize the program using the standard field-theory methods 
(first real insight into the problem was achieved by the BCS-like variational calculation ||) 
consists of introducing the field 




(7) 



P ab T AB 5 a(3 %p bl3B 

The matrices P and T defined as P = e ia A + e la S, T = (cosG + ir 3 sin Q)e"V 2 have the 
property P + P = T + T = 1. The field q operates in space of Pauli matrices abbreviated as 
rV In terms of q the Lagrangian (0) of central interest is 



3 



A) = i 

We have introduced 
where 



^-S + /i7o -A 
-70 A + 7o $—Yj — ftjo 

A = A 5 7 5 + Ay7o7 3 



q = qS (q)q. 



A s = A (1) e~ i(Q+<w cos 6 - iA (3) e~ i(<T+ ^ sin 6 



(8) 



A, 



A (4) e- i(CT+0) cos 6 - iA {2) e~ i{a+ ^ sin 6 



For rewriting any local four-quark interaction in terms of q it is convenient first to symmetrize 
all bilinear combinations with respect to if) and ip c , and then to use ip — ^(l + ^3)9, 
and^ c = ^(l-T 3 )P + T+q. 

Calculation of the quasiquark propagator requires some manual work. Abbreviating 



S(p) = 



I J 



y K L ) 

and imposing the simplifying constraint -ReA^Ay = we find 

/ = 



V p 1 /If 

7 oA + 7o(jf + + E)] 



D + 1 D 
P + M 



K 



D 



7oA + 7o(^ H 



where 



P+ = ((Po + V>),P) , D + = (Po + V) 2 ~ 4 ' 
e 2 = p 2 + £ 2 , M = [D + - (|A 5 | 2 - |Ay| 2 )]£ 



P 3 


= [D- 


h-(|A s | 2 - 


|Avf)]p 3 , 






= D + 


(Po -/*)-( 


A 5 | 2 -|Ay 


2 )(Po + /i), 


Pi 


= [D- 


h-(|A 5 | 2 + 


|Ay| 2 )]pi- 


2|A 5 ||Ay|p 2 


P2 


= [D- 


h-(|A 5 | 2 + 


|Ay| 2 )]p 2 + 


2|A s ||A v |pi 


D 


= P 2 


— M 2 = P 2 


- (P 2 + M 2 


)• 



Formulas for J and L follow from those for K and I by replacements ±/i < 
A 7oA+7 . 

For particular case Ay = considered in |§ (without E, i.e., with e p = |p |) 



(Po + /*) 



2 -e 2 



P 2 , 



) 2 + |A s l 2 l 



P 2 , 



[(A* 



) 2 + |A 5 | 2 l 



(9) 

=F/i, and 
(10) 
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We have verified (the po integration in the gap equation is trivially done using Cauchy 
theorem) that for the instanton- mediated interaction the gap equation (4.4) of || is obtained 
in an approximation of keeping only the leading term of the strongest residue. We have also 
verified that in the vacuum sector (// = 0) where the formulas greatly simplify our coupled 
gap equations for A(i) and S are in accord with corresponding formulas of [ |I~4|j . 

If the interaction is such that the solution A$ 7^ 0, X ^ is energetically favorable, the 
lowest-order self-consistent perturbation theory (nonperturbative in the coupling constants) 
turns the system of interacting massless quarks into a system of noninteracting massive 
quasiquarks with the dispersion laws given by fllPp. The corresponding specific heat exhibits 
characteristic exponential behavior. In the next mandatory step the list of the physical exci- 
tations has to be supplementented with the gapless Nambu-Goldstone collective excitations 
interacting in a calculable manner with the massive quasiquarks. 

Analogously, for the particular case As = we find the dispersion laws of the noninter- 
acting quasiquarks in the form 
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pI 



pI 



\A V ? 



4 ~ f + \ A v\ + 2 (^ - lAvf )e + (Pi + £ 2 )|A 



v\ 



1/2 



1/2 



(11) 



It nicely exhibits spontaneous breakdown of the rotational symmetry in this phase. The 
formula (|ll|) enables straightforward quantitative analysis of the Ay formation (together 
with £), apparently quite involved in the variational approach. While the dpo integration 
remains trivial, the d 3 p one becomes involved. 

Formalism presented above permits a general analysis of the phase structure of the 
deconfined QCD matter for any given interaction between quarks, i.e., not only the local 
four-fermion one. The constraint Re A $ Ay = was used merely for an easy inversion of the 
general matrix 

P-M + 2 J ReA 5 A y [p 3 7o75 - (po + A07375 + ^12] 



entering the quasiquark propagator S. 

The model is defined together with a prescription of handling the three-momentum in- 
tegrations. With the formfactor F(p 2 ) = [Is? /{p 2 + A 2 )] 1 / 2 suggested in || the dynamically 
generated quark mass falls down at large momenta within logarithmic accuracy in accor- 
dance with asymptotic freedom For obtaining realistic numerical estimates of E, A and 
T c the correct functional dependence of these quantities upon the dimensional parameters 
of the model may be crucial. 
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